Efficiency of surface-driven motion: nano-swimmers beat micro-swimmers 
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Surface interactions provide a class of mechanisms which can be employed for propulsion of micro- 
and nanometer sized particles. We investigate the related efficiency of externally and self-propelled 
swimmers. A general scaling relation is derived showing that only swimmers whose size is comparable 
to, or smaller than, the interaction range can have appreciable efficiency. An upper bound for 
efficiency at maximum power is 1/2. Numerical calculations for the case of diffusiophoresis are 
found to be in good agreement with analytical expressions for the efficiency. 

PACS numbers: 47.57.J-, 82.70.-y, 47.63.mf 



Introduction.- In recent years, much effort has been 
spent to understand biological micro-swimmers and to 
develop their artificial counterparts Mechanisms for 
propulsion of micron-scale objects include, e.g., flagel- 
lar motors, surface streaming and non-reciprocal shape 
distortions. A different class of mechanisms arc based 
on surface interactions which convert gradients in the 
environment of the swimmer into a hydrodynamic flow 
around the particle, thus propelling it forward. A 
paradigmatic example is diffusiophoretic motion in chem- 
ical gradients [U • Here, theoretical and experimental ad- 
vances have lately furnished a good understanding . 
Miscellaneous phoretic mechanisms are based, e.g., on 
gradients in temperature or electrical fields. In contrast 
to other ways of self-propulsion, they require no mechan- 
ical deformation of the swimmer. They are thus very ex- 
pedient to use for synthetic nanomotors whose mechan- 
ical degrees of freedom are hard to control. In contrast 
to the efficiency of molecular motors [l^] or other swim- 
ming mechanisms at low Reynolds numbers [ll-13|, the 
efficiency of surface-driven propulsion has, to our knowl- 
edge, received no attention so far. However, this ques- 
tion becomes relevant when energy resources are limited 
by the environment of the swimmer. Furthermore, if the 
envisioned application dictates a high swimming veloc- 
ity, the quadratic speed dependence of the dissipation 
will bring energetic aspects to attention. Industrial ap- 
plications of diffusiophoresis employ a high number 
of particles and therefore efficiency may become rather 
important here. Likewise, large dissipation could lead to 
undesired side effects like local heating. One might also 
ask from a biological perspective for the advantage of 
motility based on active processes on or near the surface 
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issues motivate us to investigate generic 
features of the efficiency of surface-driven motion, both, 
for externally and self-propelled swimmers. 

Model.- The swimmer, a spherical particle with ra- 
dius -R, is surrounded by a multicomponent fluid in a 
large container. The translational velocity of the swim- 
mer in the laboratory frame is V and it does not rotate. 
A constant, externally applied, force F m may act on the 
particle. The Reynolds number is assumed to be small 



enough that the fluid can be described by the Stokes 
equation. The variety of mechanisms which can generate 
hydrodynamic flow near the surface employ very differ- 
ent sources of energy and differ accordingly in their ther- 
modynamic description. In an isothermal steady state 
the overall energy input is given by the entropy produc- 
tion of the system multiplied by temperature plus the 
external work delivered by the swimmer. These quanti- 
ties can in principle be calculated in the framework of 
irreversible thermodynamics, thus allowing for the com- 
putation of efficiencies. Yet, the results are usually not 
analytically accessible because nonlinear field equations 
must be solved. Seeking a more general answer to the 
question of efficiency of surface-driven processes we here 
leave process-specific details aside. Instead, we concen- 
trate on the hydrodynamic efficiency eh which is a com- 
mon upper bound for the true efficiency e. The (positive) 
power output is given by P a = — VF m . The (positive) 
hydrodynamic power input Ph,i is then defined as the sum 
of power output and hydrodynamic dissipation leading to 
a bound on the true efficiency 

€ < £ h = Ti = -^7^ — TTT-^ — I 1 ) 
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where v represents the velocity and rj the viscosity of the 
fluid. E = (Vv + (Vv) T ) /2 is the deviatoric strain rate. 
We employ the customary assumption that V • v = in 
the limit of dilute solutes. This implies, besides incom- 
pressibility, that the mass densities of all fluid compo- 
nents are similar. The efficiency Eq. (JlJ depends on the 
external force F m . However, in order to find a typical 
value of eh, we focus on the hydrodynamic efficiency at 
maximum power output ej, which eliminates F m . 

Scaling of efficiency of micro-swimmers.- The 
lengthscale of the surface interactions is typically on the 
order of several nm. For particle sizes on the order of 
jUm it therefore makes sense to split the hydrodynamic 
problem into an inner problem where the fluid speed 
is strongly influenced by the surface interaction and 
an outer problem where the direct influence of the 
interaction is negligible. This is the classical boundary 
layer approximation which we will use in the following. 
We start with an estimate for power output P Q . The 
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external force in our expression for power output — F m V 
demands for the presence of a long ranged stokeslet (l7j 
and P a ~ ?7i?|V| 2 . Next, we estimate the power input 
Ph,i consisting of the hydrodynamic work rate in the 
inner and the outer region. The work rate in the outer 
region can be calculated from ordinary hydrodynamics 
and scales as r]R\V\ 2 . The dissipation in the boundary 
layer deserves a slightly more careful analysis. Here the 
thinness of the layer together with a no-slip condition 
at the surface of the particle leads to a drastic change 
of fluid velocity. This implies strong viscous dissipation. 
The thickness of the hydrodynamic boundary layer 
is denoted by L. The radial derivative is ~ 1/L and 
the speed normal to the surface is small due to the 
impermeable boundary. Hence, the dissipation rate per 
volume in the boundary layer can be approximated as 
2?yE : Vv ~ rj\V\ 2 /L 2 . When performing the volume 
integral over the local dissipation rate, the dominant 
contribution comes from a volume near the surface 
which we approximate with R 2 L. Taken together, we 
find for the power input Ph.i ~ 7y| V| 2 i? 2 / Z/ where we 
have already dropped P and the dissipation outside the 
boundary region because their relative contribution is 
0(L/R). Taken together, we find that the hydrodynamic 
efficiency scales to leading order in L/R as 

e h ~ L/R. (2) 

This generic scaling shows that any surface interaction 
whose range is considerably smaller than the size of the 
particle is inefficient in driving it. 

Upper bound on efficiency.- The above estimate for 
the hydrodynamic efficiency loses its validity for nano- 
swimmers where L is comparable to R. To derive a 
general upper limit for the hydrodynamic dissipation we 
compare the dissipation rate of the true fluid velocity field 
v with the dissipation rate in an auxiliary velocity field 
v' around a passively dragged particle. The true fluid 
velocity v is driven by any velocity independent body 
force b, arising, e.g., from surface interactions and hence 
t/V 2 v — Vp = — b. v' satisfies the homogeneous Stokes 
equation 7/V 2 v' — Vp' = and V • v' = 0. The boundary 
conditions for v' are to be the same as for v. Starting 
with the inequality (E' - E) : (E' - E) > one finds 

2? ? J E : Vvdy > J V • (v [-p'l + 2r/E']) dV = VTV 

(3) 

where T is the resistance tensor of translation, which is 
for a spherical particle given by dirnRI. Also, due to 
linearity of the Stokes equation the power output can be 
written as P Q = — V [TV + J 7 ] where J 7 is a function of 
the surface interaction forces but independent of swim- 
ming speed. Explicit expressions for T can be obtained 
[l8j |. but arc not required here. The particle velocity at 
maximum power output is V* = —T~ 1 J r /2. Using V* 



and Eq. © in Eq. flU yields 

4 < 1/2. (4) 

Remarkably, this is in a quite general sense an upper 
bound for the efficiency at maximum power of any hy- 
drodynamic motor in the Stokes regime. It is formally 
related to results for heat engines [19f but differs in the 
definition of efficiency and in that we are dealing with 
hydrodynamic systems. As demonstrated by the numeri- 
cal calculations below, the upper limit for hydrodynamic 
efficiency can be almost achieved by small swimmers. In 
the remainder of this letter we shall support these general 
considerations by a detailed treatment of diffusiophorc- 
sis where interactions with gradients of ionic or neutral 
solutes drive the particle. 

Diffusiophoresis.- As customary, the system is 
treated in the quasi-stationary limit. We employ a spher- 
ical coordinate system aligned in the 63 direction where 
r is the distance from the particle center and 9 is the 
inclination angle in the axisymmetric problem. The po- 
tential ^(r, 9) mediates interactions between the swim- 
mer and the solute concentration fields Ci (r, 9) in the di- 
lute limit. In the case of ionic solutes, symmetrically 
charged cations (ci) and anions (02) with different mobil- 
ities are present. 'J* must then be determined from Pois- 
son's equation V 2 '!' = — 47r(Ze) 2 e^ 1 [c\ — C2) where Ze 
is the charge of each ion and e is the dielectric constant of 
the fluid. The range of the ionic potential is determined 
by the Debye length I = k" 1 = [8n{Ze) 2 c°°(0)/(ekT)] 1/2 
where c°°(0]_is the concentration at r = in absence of 
a swimmer [20(. kT is the thermal energy scale. For the 
case of a non-ionic concentration gradient we use only one 
kind of solute (ci) interacting with the swimmer via an 
arbitrary, radially symmetric potential ty(r) which also 
decays on some lengthscale I. The resulting steady state 
solute fluxes are ji , j2 in the ionic, and ji in the non-ionic 
case, respectively. When neglecting convection, we have 
for solute conservation 

= V • ji, 2 = V ■ [-D 1>2 (Vci, a ± (kT)- 1 ci, 2 Vtf)] (5) 

where D\. £> 2 are the diffusion constants of the so- 
lutes. For diffusiophoresis of a passive swimmer in an 
externally maintained concentration gradient the bound- 
ary conditions are e r j-i |i- =j r = and (Vcj) |, — ^ = 
const, x S3. Also, the boundary conditions of an ionic po- 
tential are determined such that the electric current van- 
ishes at infinity (ji — J2) |r->oo = [2l|. Both, the ionic 
and non-ionic solutes mediate a body force b given by 
— (ci — c 2 )V , I , and — ciV^, respectively. Accordingly, 
the Stokes equation with ±F m |j §3 becomes 

?7V 2 v-Vp=-b, v| r=R = ! v\ r ^ 0o = -Ve 3 . (6) 

This model is only valid if the mutual interactions of 
solutes with radius a are negligible. The corresponding 



3 



corrections to the diffusion coefficient are proportionai 
to the volume fraction and we hence demand ca 3 <C 1. 
Moreover, the relative corrections of the solute-swimmer 
interactions are ~ (a/R) 2 , which should also remain <C 1 
if the solute size is to be neglected. This restricts our 
model to swimmers which are at least one or two magni- 
tudes larger than the solutes. For solutes in the A range, 
we hence require R > 30 nm. Then, the nano-swimmer 
regime corresponds to a " diffusiophoretic Debye-Hiickel 
limit" . 

Diffusiophoretic efficiency.- In order to explicitly 
confirm the scaling of micro-swimmer efficiency, Eq. ((2J , 
we employ the established theory for diffusiophoresis 
when I <C R [H, H3 • The smallness of I implies that the 
normal concentration profile near the surface is near equi- 
librium. Then ci j2 (j/, 6) sa c(9) exp [=F^(y)/fcT] where 
y = r — R and c(9) is the undisturbed concentra- 
tion of solutes outside the boundary layer. To lead- 
ing order in l/R the radial body force is compensated 
by a radial pressure change e r b « d y p which yields 
p ~ Po + kTj2i[a{y^)—c{0)]. The leading order 
contribution of the Stokes equation for lateral flow is 
rjdyVe — R~ 1 dgp « — egh. Upon insertion of the pres- 
sure and integration one obtains the boundary layer fluid 
velocity in a comoving frame Q 



ve(y,0) ~ 



kT 
V Jo 



(7) 



with the function f(y) given below |22j. Extension of the 
integral limit to infinity yields the so-called slip velocity 
v s (9) = vg(y,0)\ y ^oo at the interface between bound- 
ary layer and the outer flow. To calculate the particle 
speed V from the slip velocity one matches v s as bound- 
ary condition to an outer hydrodynamic solution where 
b = 0. The result is V « ^ R + \ /* v 3 {9) sin 2 6d6. 
The efficiency is determined by inserting the inner solu- 
tion Eq. (J7J into Eq. ([T]) while the matched outer solution 
does not contribute in leading order. The result for the 
efficiency at maximum power becomes 



6h ~ 32f*{d e c{9)) 2 sm9d9 R 
where we have defined the dissipation length 

r U °°yf(y)M 2 



^J oa f(y)^y'f(y')dy'dy- 



(8) 



(9) 



The dissipation length Ld is a measure for the radial ex- 
tension of the layer where velocity gradients are strong. 
It is expected to be similar to the thickness of the layer 
where the fluid is driven by the body forces. For ionic 
solutes Ld is to leading order proportional to the Dcbyc 
length k" 1 [H. 

In the case of non-ionic solutes the expression for Ld de- 
pends on the choice of 4 r (r). Here it is of interest to com- 
pare Ld with the thickness of the layer of excess solute 



L* 0, [24[ because the latter can be inferred indirectly 
by measuring the diffusiophoretic speed. It is usually on 
the order of 10 nm We find L D = L* = 1/2 in 

the concrete cases of a hard-core repulsion as well as for 
^{y)/kT = {l/y) n with n — > oo and small I. Moreover, if 
one replaces the inner integral limit y in the denominator 
of Eq. (|9]) with oo then Ld = L*/3. This, together with 
evaluations of Eq. ([9]) for various ty(r) shows that it is 
safe to estimate the magnitude of Ld through |L^| in the 
non-ionic case. 

We have also investigated the effect of hydrodynamic 
slip boundary conditions [26j on e£. These reduce the 
hydrodynamic dissipation in the boundary layer but ad- 
ditional dissipation occurs directly at the surface. Effec- 
tively, Ld is increased by a few nm. A good absolute 
value for efficiency can, however, not be achieved in this 
way when Ld -C R is still valid. 

Numerical analysis.- In order to extend our analysis 
to nanoparticles where I > R we solve the coupled dif- 
fusion and hydrodynamic equations numerically In the 
case of ionic solutes the nonlinearities are avoided by ex- 
panding the solution for low dimcnsionless surface charge 
density q = AirZeQ^s k/cT) -1 as done in [20| . The con- 
stant q can then be related to the electrostatic potential 
on the surface Hf{R) = kT q kR (1 + kR)^ 1 +0(q 3 ). For 
non-ionic solutes we can directly solve Eqs. ([5]) and ((6]) 
after choosing a potential Vf'(r). 

The resulting efficiency at maximum power e£ is dis- 
played in Fig. [TJ Note that e£ is independent of 77 
and (Vcj) |r->oo- The absolute concentration level en- 
ters only for ionic solutes via the Debye length. For 
Ld < i? we find good agreement of the numerical re- 
sults with the boundary layer theory where Eq. © pre- 
dicts ejj = (3/8) Ld since the concentrations C\${r, 9) are 
linear in cos^. For Ld > R the efficiencies of nano- 
swimmers are actually in the range of the upper bound 
given by 1/2. This is due to the low dissipation when the 
swimmer moves in a spatially slowly varying flow field. 
The deviations from the analytical predictions are most 
pronounced in the intermediate regime of Ld — R- For 
ionic solutes (inset of Fig. Q]) we find that the efficency 
is larger when the two solutes have different mobilities 
v 7^ 0; see (23- Then electrophoresis in self-generated 
concentration disturbances around the particle influences 
the swimming. 

Janus particles.- To complement the results for pas- 
sive swimmers we also investigate active Janus particles 
where the neutral solute concentration gradient is main- 
tained by a chemical reaction on the surface 0-0]. This 
is modeled through the boundary conditions of Eq. ([5]) 
by e r ji| r= fl = a(l + cos9) and ci\ r ^oo = const, where 
a is the effective solute production rate per area. The ef- 
ficiency £j* becomes independent of a and differs only in 
the regime of Ld — R from the results for swimmers in an 
externally applied gradient. Active Janus particles are in 
this regime hydrodynamically less efficient (see Fig. [1} . 
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This relates to the fact that long range potentials are not 
so effective in driving the fluid when the concentration 
gradient decays to zero away from the particle. A prob- 
lem occurring for small Janus particles is the fast con- 
centration field homogenization through their rotational 
diffusion on the timescale of 4ttt]R 3 (fcT) -1 . It may thus 
be necessary to fix such motors directionally. 
We also mention that the full efficiency e of Janus parti- 
cles has to take into account losses due to maintenance 
of spatial concentration gradients and due to chemical 
reactions. The power input then reads Pi = P^h — 
/ ji'V/x» + J2k Ak Tk) dV where ^ is a chemical po- 
tential of species i including Vfj. and are affinity 
and rate of the k-th. chemical reaction. 



been explored. In this letter we make an attempt in this 
direction by suggesting a generic scaling relation for the 
efficiency of surface-driven motion. It provides a widely 
applicable and simple concept to estimate hydrodynamic 
emcency without detailed knowledge of the system. Fur- 
ther, we show with analytical and numerical calculations 
that phoretic nano-swimmers offer energetic advantages; 
in particular with ionic solutes, where the Debye length 
can be tuned. Taken together, we see inspiring perspec- 
tives for artificial nanomotors which, reminiscent of ac- 
tual biological motors, could possibly move not only in a 
controllable but also in an efficient way. 
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FIG. 1. Efficiency at maximum power e£ for diffusiophore- 
sis in non- ionic solutes. (□): Van der Waals attraction 

,3 r>3„,-3/ 



9(y) = -W/9Aa i R i y- <s (y + 2R)~ 6 where A = lfcT and 
a is the solute radius. a/R — [l(F 4 ...0.l] and attraction 
is truncated at y = a. Concentration gradient is established 
externally. (■): Same ^(y) as before but concentration gradi- 
ent is produced by a Janus particle. ( ): Generic repul- 
sion *(y) = kT exp {-y/l) with l/R = [lCT 4 . . . 10 2 ] . Con- 
centration gradient is established externally. ( ): Same 

ty(y) as before but concentration gradient is produced by a 
Janus particle. Inset: difhisiophoresis in an externally estab- 
lished gradient of ionic solutes and n' 1 /R = [l(T 4 . . . fO 2 ] . 

(• • •): q = 0.1, v = 0. (o o o): q = -0.5, v = 10. ( ): 

q = 0.1, v = 1. 

Swimming speed. - The energetic differences of micro- 
and nano-swimmers are accompanied by differences in 
the swimming speed. As a simple example, for micro- 
swimmers in neutral solute with F m = one has V 

f ^ 2 (v Cl )| r 



0. However for swimmers with I 3> R 



we find the scaling V ~ ^ ^ (Vci) | r ^oo [27J. Hence, 
aside from efficiency, nano-swimmers are qualitatively 
different from micro-swimmers in that their size can mat- 
ter for their mobility. 

Conclusion.- Although phoretic effects are known for 
more than a century their energetic aspects have hardly 
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